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Introduction 

The  purpose  of  this  technical  report  is  tuuof old: 

f  to  critically  examine  the  ITDfl  FiSlU  barrier  and  area  search 
models  for  mathematical  accuracy  and  modelling  reasonableness,  and 
2.  to  report  on  the  use  of  ITDfl  programs  on 
USS  CARL  VINSON  (CVN  701  during  a  id  day  period  of  high  intensity, 
exercise  operations  (RIMPflC  86). 

Available  for  examination  mere  ITDfl  Ver.  1  (including  the  BASIC 
language  source  code),  the  more  recent  ITDfl  Ver.  2  (without  source 
code),  and  the  JOTS  II  Technical  Reference  Manual  (Ref.  I).  Reference  is 
made  in  this  report  to  Version  1  source  code  to  explain  the  functioning  of 
the  Version  2  program.  This  mas  required  because  the  Version  2  source 
code  mas  unavailable.  Fortunately,  the  barrier  and  area  search  modules 
appear  functionally  very  similar,  if  not  identical,  in  both  ITDfl  versions. 

So  the  lack  of  Version  2  source  code  is  not  judged  to  be  especially 
critical  to  this  examination. 


1 


fiSLU  Barrier  Search  Model 


Brief  Description  of  Model.  The  ITDR  RSUJ  barrier  search  model  is 
accurately  described  in  Fief.  1.  Its  description  here  is  primarily  to 
establish  notation. 

The  barrier  search  model  estimates  the  probability  of  detection  of  a 
target  submarine  penetrating  a  back-and-forth  RSUU  barrier  (Figure  1.). 
The  barrier  penetration  is  considered  in  "target-stationary  relative 
space".  That  is,  the  target  is  assumed  to  be  stationary,  and  all  relative 
speed  for  the  encounter  is  provided  by  the  searcher.  The  searcher's 
relative  speed  component  across  the  barrier  front  is 

Vs  costa), 

where  Vs  is  searcher  speed,  and  a  is  the  barrier  advancing  angle  (when 
,.v>0]  or  retiring  angle  (when  a<0).  The  component  of  relative  speed 
perpendicular  to  the  barrier  is 

V  T  +  V§  sin(cd, 

where  VT  is  target  speed. 


I 

Target  window,  UJ 
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Figure  1.  Barrier  Geometry 


For  barrier  penetration  points  uniformly  distributed  across  a  target 
window  (UJ),  the  probability  of  the  target  coming  within  range  R  of  the 
searcher  (PD)  is  the  ratio  of  the  area  “covered"  by  the  searcher  in 
relative  space  to  the  total  area  that  could  be  occupied  by  the  target.  In 
Figure  2.  the  covered  area  is  shaded.  In  the  ITDfi  barrier  search  model, 
this  area  is  obtained  by  numerical  integration.  Allowance  is  made  for 
sensors  which  sweep  less  than  360°  by  reducing  the  detection  range  R. 
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T  T 
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Figure  2.  Calculation  of  P0 


Discussion.  Although  the  modelling  concepts  used  are  good.,  there  are 
some  minor  problems  in  the  implementation. 

I.  The  calculated  F'D  is  too  large  for  target  luindotu  sizes  (UJ)  less  than 
2(B/2  -  RK).  This  problem  results  because  the  limits  of  the  numerical- 
integration  are,  in  some  instances,  improperly  set.  Specifically,  in  line 
51970  of  the  ITDA  Ver.  1  BASIC  code,  the  upper  limit  of  integration  is 
B/2  -  RK.  It  should  be  min  (B/2  -  RK,  IJJ/2).  With  the  original  limits,  the 
integration  can  be  performed  outside  the  target  window,  W.  This 
increases  the  calculated  F'D  above  the  proper  value.  To  illustrate  with 
an  example,  let 
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R  =  detection  r  ange  =  10  nrn 
Vs  =  searcher  speed  =  19  knots 
V'T  =  target  speed  =  19  knots 
B  =  barrier  midth  =  89  rim 
IJJ  =  target  window  =  19  nrn 
oc  =  advancing/retiring  angle  =  9 

K  =  kinematic  enhancement  factor  =  (1  +  I. V 3/ V j j ^  1 1  ~  =  1.414 
The  ITDfi  model  gives  a  probability  of  detection  of  I.  However,  for  any 
UJdB  -  2RK),  the  correct  answer  is  2RK/B  -  .35. 

2.  PD  is  improperly  calculated  when  either  Vs  cos(ct)  or  VT  -  Vs  sin(a) 
is  negative.  The  cause  of  this  problem  is  that  the  integrand  in  the 
numerical  integration  in  line  51169  of  the  BASIC  code  is  allowed  to 
become  negative  whenever 

Ksqr  =  Vs  costa)  /  (VT  +  Vs  sinta))  <  9. 

Adding  negative  terms  in  the  numerical  integration  results  in  the 
calculated  PD  being  too  small.  To  illustrate,  continue  the  previous 
example  with  IU  =  199  nrn  and  a  =  -19°.  Then  PD  is  correctly  calculated  as 
.31.  Setting  a  =  - 179°  should  give  the  same  result,  but  yields  .29. 

The  failure  becomes  particularly  spectacular  when  the  denominator 
of  Ksqr  is  exactly  9.  For  example,  setting 
R  =  19  nrn 
V's  =  29  knots 
V t  =  19  knots 
B  =  89  nrn 
ill  =  199  nrn 
a = -39° 

results  in  a  PD  of  1.42  x  I9b 
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R  simple  fix  for  this  problem  is  only  allow  <>.  such  that 
~90°  ■■  ct  90°  arid  Vj  +  vg  sin f oc }  s  0. 

The  other  cases  are  either  redundant  or  result  in  the  searcher  moving 
"south"  faster  than  the  target. 

3.  Sweeping  sensors  are  modelled  in  an  unusual  manner  The  user  is 
given  the  option  of  calculating  P0  under  the  assumption  that  the 
detecting  sensor  sweeps  ±8°  either  side  of  a  center  bearing.  The 
effective  detection  range  then  becomes  Rsin(minl6,90°)).  The  problem 
with  this  approach  is  that  it  assumes  that  the  center  bearing  is  along  the 
relative  track  of  the  searcher.  Rn  more  reasonable  assumption  is  that 
the  center  bearing  be  along  the  searcher's  absolute  (i.e...  geographic) 
track.  See  Figure  3. 


Figure  3.  Orientation  of  Swept  Hrea 


These  two  assumptions  lead  to  significantly  different  calculated 
probabilities  of  detection.  For  example,  in  the  limit  as  8  and  ct  decrease 
to  0,  the  ITDR  model  gives  a  probability  of  detection  of  0.  Whereas 
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assuming  that  the  sensor  is  looking  directly  ahead  along  a  fixed  line  of 
bearing  gives  a  probability  of  detection  of  mini  I,  Fl/UJ).  Note  that  in  this 
case,  there  is  no  increase  in  Pq  as  the  searcher  increases  speed  li  e  ,  no 
kinematic  enhancement).  The  probability  of  detection  is  simply  that 
fraction  of  the  possible  barrier  penetration  points  that  are  covered  by 
the  finite  detection  "ray"  of  the  sensor.  For  nonzero  values  of  8  and  ot, 
similar  closed  form  expressions  can  probably  be  derived,  although  the 
geometry  becomes  more  complex  and  performing  the  numerical 
integration  may  be  preferred. 

4.  This  is  a  comment,  not  a  problem.  For  c*=0°  and  a  369°  sensor,  it  is 
possible  to  approximate  the  ratio  of  the  shaded  area  to  the  total  area  in 
Figure  2.  without  numerical  integration.  True,  the  HP9929  is  not  being 
computationally  stressed  by  the  integration,  but  there  is  something  to 
be  said  for  the  simplicity  of  closed  form  solutions.  Using  the  notation 
introduced,  the  probability  of  detection  given  a  barrier  penetration  is 
approximately 


min  (2RK/B,  1J  for  LU  $  B  —  2RK 


Po  = 


1  -  (1/UJB)(max  (,5(B  +  UJJ  -  RK,  0jj2  for  B-2RK  <  UJ  <  B  +  2R 

((B  +  2RJ/UJ)  (1  -  (max  (B  -  RlK-1),  0)J2  /(B2+2RBJ)  for  IJJ  >  B  +  2R 


where  K  =  (I  +  (V3/VT)2),/2  as  before.  To  obtain  this  expression  for  F'0, 
the  circular  boundary  of  the  shaded  area  in  Figure  2.  is  replaced  with  a 
straight  line  The  resulting  region  has  an  area  which  is  greater  than  or 
equal  to  the  area  of  the  original  shaded  region  So  this  calculation  gives 
an  upper  bound  for  P0.  Figure  4.  is  a  comparison  plot  of  PD  calculated  by 
the  ITDfi  model  and  the  approximation  above  The  ITDfi  model  values 


greater  than  .35  are  in  error  and  result  from  improper  limits  of 
integration  used  when  calculating  PD  for  small  IU. 


Figure  4.  Comparison  of  ITDR  Model  Results  and  Approximation 


5.  This  is  a  second  comment.  One  advantage  of  using  numerical 
integration  for  evaluating  PD  is  that  it  can  accommodate  a  nonuniform 
distribution  of  target  penetration  points.  This  advantage  is  not 
exploited  in  the  current  barrier  search  program.  For  example,  if  the 
target  knows  that  the  barrier  is  in  place,  then  an  attempt  might  be  made 
to  penetrate  near  the  one  of  the  ends  of  the  barrier  to  minimize  F'0.  If 
numerical  integration  is  used  to  evaluate  PD,  then  it  is  a  simple  matter  to 
weight  the  possible  barrier  penetration  points  according  to  any  specified 
probability  distribution.  Specifically,  if  the  barrier  penetration  point  x 
has  a  probability  density  function  g(x)  and  PD(xl  is  the  probability  of 
detection  given  a  barrier  penetration  at  x,  then 


U 


pd  =  J‘  pd(x)  gtxi dx 

HI  I  ;< 

s.t  q-'x»0 

6  R  third  comment.  R  possible  problem  ujith  this  advancing/retiring 
barrier  methodology  is  that  the  detection  probabilities  are  conditional  on 
the  target's  penetration  of  the  barrier  Changing  the  angle  of  searcher's 
advance  or  retreat  ijuill  affect  both  the  probability  of  detection  given 
barrier  penetration  and  the  probability  of  barrier  penetration.  The 
second  effect  is  not  accounted  for  in  this  model.  For  example,  a 
retreating  barrier  can  dramatically  increase  the  probability  of  detection 
given  a  barrier  penetration,  but  it  also  reduces  the  probability  of  such  a 
penetration  during  any  specified  time  interval.  In  the  limit  as  the 
searcher  speed  perpendicular  to  the  barrier  increases  to  the  target 
speed,  the  predicted  PD  increases  to  1.  However,  the  target  can  not  now 
overtake  the  retreating  barrier,  so  the  probability  of  a  barrier 
penetration  occurring  at  all  becomes  zero.  This  line  of  argument 
suggests  that  the  calculated  probabilities  of  detection  may  be  too  large 
for  retiring  barriers  and  too  small  for  advancing  barriers. 

flSUJ  Area  Search  Model 


Brief  Description  of  Model,  fls  described  in  Ref.  1,  probability  of 
detection  for  an  area  search  is  given  by 
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PD  =  I  -  expl-ZTi  VRt/A), 

inhere 

R  =  detection  range 
fl  =  search  area  size 
t  =  time  alloiued  for  search 
V's  =  searcher  speed 
V T  =  target  speed, and 
VR  =  mean  relative  speed 

2p 

=  11/2  rr  I  j  (V.-/  +  VT2  +  2VsVtcos  4>)1/2  d<t> 
u 


(2 


Discussion 


1.  The  mean  relative  speed  calculation  can  be  made  more  efficient. 
In  the  original  BASIC  code  and  in  Ref.  1,  the  calculation  of  VR  is  as  shown 
in  Equation  (2).  A  preferred  expression  for  numerical  evaluation  is 

Tl 


(1/  JT  )  | 


-  2VsVtcos  4') 1/2  d$  . 


(3) 


Equation  (3)  makes  use  of  the  symmetry  of  the  cosine  function  to  reduce 
the  extent  of  the  numerical  integration,  and  has  the  correct  sign  for  the 
2VsVt  cos'4  term  in  the  Lam  of  Cosines.  These  problems  are  to  some 
extent  cosmetic,  as  cosine  symmetry  causes  both  (2)  and  (31  to  give  the 
same  numerical  results. 

It  is  noted  that  VR  can  be  approximated  mith  the  following  closed 
form  expression: 

VR/VS  ~  .36  max(VT/Vs  ,  1)  +  .64  (1  +  (VT/'VS)2)1/2  (4) 

Figure  5  contains  plots  of  the  estimated  VR/VS  versus  VT/VS  and  a 
multiplicative  correction  factor  to  convert  the  estimated  values  to  true 
values.  The  estimate  is  always  within  1.1%  of  the  true  value. 
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SPEED  RATIO  (TGT.  SPD./SEARCHER  SPD.) 


Figure  5  Estimated  Relative  Speed  Ratio  vs  Speed  Ratio 


2.  R  Monte  Carlo  area  search  simulation  mas  conducted  to  test  the 
accuracy  of  Equation  11).  The  search  geometry  is  shown  in  Figure  6  The 
searcher  starts  in  the  lower  Left  corner  of  a  IRQ  rirn  x  IQ 9  nrn  search  area 
arid  conducts  a  systematic  search  at  speed  Vs  with  a  cookie-cutter 
sensor  having  a  range  of  IQ  nrn.  The  target's  starting  position  is 
uniformly  distributed  over  the  search  area.  The  target  conducts  a 
random  tour  at  speed  VT  with  a  mean  time  between  course  changes  of  1 
hour.  That  is,  the  target  track  consists  of  connected  Line  segments, 
each  of  which  has  a  course  selected  from  an  independent,  uniform 
distribution  between  0°  and  360°.  find  the  time  on  each  leg  is  selected 
from  an  independent,  exponential  distribution  with  mean  1  hour.  When 
the  target  encounters  a  search  area  boundary  a  perfect  reflection  is 
performed. 


1 00  nm 


Figure  6.  Simulation  Geometry 


Seven  simulations  mere  conducted,  each  mith  Vs  and  VT  selected  so  that 
the  mean  relative  speed,  as  calcuated  from  Equation  13),  mas  15  knots. 
The  results  are  plotted  in  Figure  7,  together  mith  P0  calculated  from 
Equation  I'D.  The  fit  betmeen  (1)  and  the  simulation  mas  reasonable  when 
Vc  and  V T  mere  nearly  equal.  Otherwise  the  fit  was  not  especially 
impressive. 

The  simulation  results  were  sensitive  to  the  mean  time  between 
target  course  changes  and  the  searcher  track  used  (neither  of  which  are 
represented  in  Equation  (!)].  In  particular,  the  simulation  showed  that  a 
search  track  starting  near  the  center  of  the  search  area  gives  better 
search  performance,  especially  for  small  values  of  Vs. 

3.  The  need  for  an  area  search  simulation.  For  this  area  search 
example.  Equation  (D  does  not  model  a  systematic  search  for  a  randomly 
moving  target  very  accurately,  find  unfortunately  the  search  literature 
has  little  else  to  offer  in  the  way  of  closed  form  expressions  for  PD 
when  the  target  track  is  random  and  the  searcher  conducts  a  systematic 
search  (e.g.,  parallel  sweeps,  spiral  in,  or  spiral  out).  This  suggests  that 
some  sort  of  numerical  solution,  probably  a  Monte  Carlo  simulation,  may 
be  required  to  assess  the  effectiveness  of  such  an  area  search.  This 
simulation  could  be  very  general.  It  could  allow,  for  example: 

1.  different  areas  for  the  target  and  searcher, 

2.  a  customized  target  motion  model  (e.g.,  patrol  or  transit) 

3.  any  desired  search  track, 

4.  multiple  searchers  and  targets, 

5.  counterdetection  by  the  target. 


PROBABILITY  OF  DETECTON 


Mean  Relative  Speed  is  1  5kts 
Mean  Rate  of  Course  Changes  is  1/hr 


o 


Figure  7.  Probability  of  Detection  by  Simulation 
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Output  could  include  the  searcher  trad  Is)  and  graph  of  probability  of 
detection  by  time  t.  With  a  cursor  or  light  pen  input  of  area  boundaries 
and  track  points,  problem  setup  time  cuould  be  reduced  This  simulation 
might  be  of  particular  use  for  battle  group  flSW  planning,  since 
frequently  several  FiSlU  platforms  search  for  the  same  target 
submarine 


fit  Sea  Observations 


Background  During  the  period  2-18  June  1986,  the  author,  as  a 
commander  in  the  Naval  Reserve,  mas  attached  to  the  staff  of 
Commander  Carrier  Group  Three,  fit  that  time,  the  staff  mas  embarked 
on  USS  CARL  VINSON  (CVN  70)  and  operating  in  support  of  the  exercise 
RIMPRC  86. 

RIMPfiC  86  mas  a  two-carrier  battle  group,  multinational  exercise 
involving  approximately  60  ships  and  30,000  at-sea  personnel  The 
purpose  mas  to  demonstrate  the  fleet's  ability  to  operate  effectively 
against  a  realistic  air,  surface,  and  subsurface  threat.  During  this 
exercise,  use  of  the  ITDH  system  by  the  carrier  group  staff  and  CARL 
VINSON  ship's  company  mas  observed 

Exercise  Scenario.  The  general  scenario  mas  that  a  coup  had  toppled 
the  pro-West  government  of  White  Country  (represented  for  purposes 
of  the  exercise  by  the  Hawaiian  Islands)  The  US  naval  response  mas  to 
send  two  carrier  battle  groups  (Ranger  and  Carl  Vinson)  to  White 
Country  as  a  shorn  of  force.  These  battle  groups  mere  opposed  by 
some  60  ships  simulating  Soviet  platforms.  As  the  exercise  progressed. 
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the  battle  groups  approached  iJJhite  Country  close  enough  to  launch  air 
strikes.  Opposing  surface  action  groups  and  submarines  engaged  as  the 
battle  groups  came  into  range. 

JOTS  and  DOTS,  The  ITDR  system  on  CRRL  VINSON  is  called  VINSON 
DOTS,  which  appears  to  be  a  one-of-a-kind  installation,  funded  by 
ORRF'fi,  through  NOSC,  San  Diego.  There  was  no  formal  documentation 
for  this  system  aboard  CRRL  VINSON,  but  the  overall  program  structure 
and  functional  modules  seem  very  similar  to  the  older  JOTS.  In  fact, 
much  of  the  BRSIC  code  has  been  carried  over  directly  from  JOTS  (as 
indicated  by  documentation  comments  in  the  BRSIC  code).  The  principle 
DOTS  enhancements  included  the  organization  of  the  contact  database 
and  the  implementation  of  automated  input  of  LINK  11  data. 

The  controlling  terminal  for  VINSON  DOTS  is  an  Hewlett-Packard  9020 
located  in  the  UJar  Room  of  the  Flag  Spaces.  Feeding  this  station  is 
another  HP  9020  in  the  Carrier  Intelligence  Center  (CVIC).  This  unit 
processes  and  correlates  intelligence  information  generated  off  the 
ship.  The  correlator-tracker  used  is  POST  (Prototype  Ocean  Surveillance 
Terminal).  Also  feeding  into  the  Flag  HP  9020  is  real-time  NTDS  (Link  11) 
data.  The  Flag  terminal  displays  and  operates  on  available  contact 
data  and  also  sends  this  data  to  HP  9020  machines  in  the  RSLU  Module 
(RSLUMQO)  and  Combat  Direction  Center  (CDC).  Another  HP  9020  is  located 
in  the  Tactical  Operations  Plot  (TOP)  which  runs  JOTS  III  and  is  currently 
not  in  the  VINSON  DOTS  network  (and  so  does  not  have  access  to  real¬ 
time  NTDS  information). 


The  VINSON  DOTS  network  accepts  data  only  from  the  Flag  terminal 
and  these  data  currently  come  from  three  sources:  II  manual  input,  2J 
NTOS  input,  and  3)  POST  input. 

The  NTDS  input  is  automated  and  requires  no  human  filtering 
Automated  input  is  also  available  from  POST,  but  due  to  the 
classification  of  some  POST  data,  this  feature  is  not  presently 
implemented.  To  remove  the  more  sensitive  data,  POST  contact 
information  must  be  manually  sanitized  and  then  entered  into  the 
VINSON  DOTS  database  from  the  Flag  terminal.  Automatic  sanitization 
and  entry  of  POST  data  is  a  planned  enhancement. 

Even  though  contact  data  to  be  shared  over  the  network  must  be 
entered  from  the  Flag  9Q20,  any  other  terminal  may  be  used  to  enter  or 
modify  contact  information  to  be  used  exclusively  on  that  terminal. 

This  allows  users  to  customize  their  databases. 

Use  of  DOTS  Aboard  CVN  78.  Of  the  many  tactical  decision  aids  that  DOTS 
provides,  the  principle  observed  uses  during  RIMPAC  86  were. 

1.  Maintaining  an  contact  history  database  for  the  exercise, 
containing  all  surface,  air,  and  subsurface  contacts. 

For  each  detection  recorded,  some  or  all  the  following  data  were 
saved:  time,  location,  target  course  and  speed,  detecting  sensor  and 
detection  parameters  (dimensions  and  orientation  of  elliptical  spas  or 
bearing  spread  and  inclusive  ranges  for  bearing  boxes).  In  addition, 
■JOTS  III  (but  not  DOTS)  provides  a  comment  field  which  was  used 
extensively  by  the  TOP  watchstanders  (responsible  for  maintaining  the 
surface  plot). 

2.  Plotting  current  or  past  contact  positions. 


There  were  two  principle  ways  this  capability  was  used: 

a.  Displaying  the  current  position  of  some  group  of  contacts  le  g.. 
Blue  surface  ships),  or 

b  Displaying  the  contact  history  (i.e.,  track)  of  a  few  contacts. 

Color  and  NTDS  symbology  were  effectively  used  to  identify  the 
platform  type  and  threat  status  (hostile,  friendly,  or  unknown). 

The  processing  and  display  of  Link  11  data  transfer  was  fast  enough 
to  follow  the  air  war.  This  was  demonstrated  during  several  air  raids 
on  the  battle  group  where  actual  aircraft  were  used,  (fl-3's  simulated 
bombers  with  A-4's  and  fi-7's  flying  as  air-launched  cruise  missiles.) 


-  e.g.,  if  BATES  was  underway  from  San  Diego  9890  on  2  June,  could 
she  be  in  the  Hawaiian  Operation  area  by  9809  on  10  June  with  a  15  knot 
speed  of  advance? 

r  with  a  18  knot  speed  of  advance  going  through  specified  points, 
where  will  the  contact  be  in  24  hours? 

4.  Plotting  spas,  patrol,  areas,  4-lUhiskey  grid  (i.e.,  screen) 
assignments 

5.  Preparing  "executive  summaries"  of  the  contact  situation. 


No  direct  use  of  the  ARID,  ASUI  or  ELU  decision  aids  was  observed, 
although  several  of  the  Carrier  Group  Three  staff  members  were  very 
knowledgeable  in  their  use. 


0  b  s  e  i '  v  9  d  P  ro  b  1 e  rn  s . 

1.  Resolution  on  the  HP  9828  screen  was  insufficient  to  display  many 
contacts.  (The  POST  display,  using  an  HP  98280  machine,  was  much 
better.) 

2.  TOP  watchstanders  used  JOTS  III  rather  than  DOTS  because  JOTS 
III  is  more  complete. 

3.  The  HP  Think  Jet  printer  generally  only  provided  marginally 
acceptable  hard  copy.  The  installed  thermal  printer  on  some  of  the 
HP  9828  machines  provided  a  finer  resolution  output. 
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